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Abstract 

In recent years, researchers have discovered various excited states of electrons having an excited 
domain with a structural order different from the original one. We examine a dynamics of elec- 
tron system and phonon system by using the time-dependent Hartree-Fock approximation. We 
examine time evolution with spin degree of freedom. We found a thermally excited paraelectric- 
to-ferroelectric crossover. The metastable state is spin pierls state. We also examine the Landau 
expansion of electron system at quarter-filled. 
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The paraelectric-to-ferroelectric crossover caused by electrons displacement is induced 
by an electron in an paraelectric insulating solid when it is excited by a electric external 
field.- The original concept of this relaxation assumes this crossover to be a microscopic phe- 
nomenon. The photo induced paraelectric-to-ferroelectric crossover can be called a photo- 
induced structural phase transition. In this paper, we examine the theoretical studies on 
photo induced problem. Photo-induced phase transitions have attracted much attention. 
This structural phase transition has often been described as a domino effect. It is highly 
nonlinear and cooperatives, in particular its time evolution is many body problem. The 
photo-irradiation a large energy is imported by the photo induced irradiation that is not 
accessible in thermal equilibrium. Thus even if photo-induced phase is one that is realized 
in such conventional method, the mechanism of the transition is not trivial at all. A large 
enhancement of dielectric constant by UV light irradiation is reported for the first time in 
quantum paraelectric family of perovskite oxides, such as SrTiOa and KTiOs, under a weak 
DC field.- The photo irradiation induce the paraelectric-to-ferroelectric crossover, as a func- 
tion of temperature. Without photo irradiation, the dielectric property of these material 
is Quantum paraelectric. The DC electric field suppress the quantum fiuctuation, which 
enhance the quantum domino catastrophes. 

Photo-induced transitions have been studied in many method. Whether a deterministic 
or stochastic approach accurately describes the time evolution of electron system would 
depend on the time scale of observation. The recent study obviously needs to be described 
in a deterministic manner.^ Their charge and lattice motion is calculated by using the time- 
dependent Schrodinger equation. In this paper, we study another itinerant electron system 
and phonon system. The present study differs from other studies, because we treat the 
electron transportation and many body interaction. 

We examined the time evolution from the initial state to a stationary state of the electron 
system with random initial condition. During the evolution, the excited state changes into a 
metastable state. This method has advantage because in Ref."^ , the spin degree of freedom 
can not be treated. Though, in our method, we can not treat the eigenstate in stead of that 
we can examine time evolution with spin degree of freedom. We also examine the Landau 
expansion of electron system at quater-filled. 

The electron system and classical phonon system is described by a time dependent 
Hartree-Fock approximation. The classical phonon system is evaluated by leap-fiog 
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method. We assume a that Hubbard model for thermal-fluctuation-induced paraelectric- 
to-ferroelectric crossover, as follows, 



Hei = Ei,^[/if(4^Q+i,^ + /i.c.) (1) 
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where to is tight binding hopping strength, and U is on-site Coulomb repulsion, cj^^, Ci^a 
is creation and annihilation operator at cite i, spin a. Hp^ phonon system Hamiltonian, 
Pi is momentum, M is mass k is stiffness, u[i] is distortion and a is dimerization coupling 
constant between electron system and phonon system. We will examine the simple Hubbard 
model. While, the complexity of degrees of freedom are different, the universality should 
be the same, because the important quantity is the single band electron system. The time 
evolution is described by the correlation function, 

g[z,j,tr' =<cl{t)cj,At)> ■ (4) 
The time evolution equations of a spin diagonal correlation function are 
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However, there is no eigcnstatc for this correlation function, in stead we can examine the 
time trace of system with random initial condition numerically. The equation corresponds 
to the random phase approximation. In present study, d-dimensional quantum field theory 
is equivalent to 2d classical field theory. During the time evolution, the initial excited state 
changes into the local ground state. 

The Bloch- Winger representation of the correlation function is, 

g'[x,p]^J2^'"''9''[k,l], (6) 

m 

where x — k — I, m — {k + l)/2. This representation is not positive definite, therefore 
probabilistic interpretation can not be applied. 



First, we examine numerical solution. The numerical method comprises the Crank- 
Nicholson method^ and the Gauss-Seigel Method. The time evolution of the Schrodinger 
equation is described by the unitary condition to 2nd-order approximation: 

1 + i^Ht 

exp(-iHt) ~ 1 . (7) 
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In present numerical study, we set the number of electrons to be 4 due to quarter-filled 
system. Let us examine the time dependence of the electron order parameter. CDW is 
the intensity of the alternative derivative of local charge density. AF is the intensity of the 
alternative derivative of local spin density. The order parameters are defined by 

- 2(7[^,2,r), 

(8) 

For the whole system the physical parameter is as follows, to = 0.1, f/ = 0.5, = 8, M = 500, 
a = 0.1 and unit of the time is normalized to the above parameters. The initial conditions 
are (^[i, j, t]'^ = 0.25 + /? * n where n is complex random value with absolute value 1 for all 
i,j.cr. and are jSei = 0.01 or jSei = 0.1. This means the initial condition is a quarter-filled 
homogeneous state with off-diagonal coherence. For the Bloch-Wigner representation, after 
the time evolution, the momentum distribution of the charges is oscillating and has domain 
in real space. The strengths of higher moment are weak. This result is similar to the Haldane 
conjectures. 

The CDW state order parameter is 4k p CDW order parameter.- The AF state order 
parameter is Akp ant ferromagnetic order parameter. Thus, the Akp CDW is quantum 
paraelectric order parameter. The dimerization enhance the Akp CDW state. The AF state 
order parameter represent spin pierls instability. For strong noise, AF state order parameter 
decreases and CDW state order parameter increases with time. The time trace of the two 
order parameters become in phase. The time trace of order parameters shows the phase of 
two order parameters depends on noise strengths. Exactly the final stage, CDW state order 
parameter is positive for the figs.l and 2, and AF state order parameter phase is become 
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FIG. 1: Time dependence of electron order parameter. /?ez = O.Oland 0.1 no phonon 

negative for figs. 3 and 4, when strong randamness. This phase locking due to the quasi-one- 
dimensional interchain coupling is assumed. The figs.l and 2 is ground state and figs. 3 and 
4 is metastable state. The strong randomness leads to metastable state. 

Due to complex initial condition, the transition from metastable state to the order ground 
state occurs. When /5=0.01, the AF order parameter increases, this indicates the Wigner 
crystal insulator. When /5=0.1, the AF order parameter become negative, the cross over to 
spin pierls state (ferroelectric) occurs. An experiment has shown the existence of a photo- 
induced paraelectric-to-ferroelectric phase crossover.- Thus, we have examine this crossover 
from the viewpoint of non-equilibrium condensed matter theory by using a Hubbard model 
Hamiltonian in view point of time dependent Hartree-Fock calculation. The transverse photo 
effect behaves as thermal fluctuation. The thermal fluctuation creates a non-equilibrium 
state with finite order parameter and settled random initial condition. The complex random 
initial condition correspond to the excited phonon with imaginary self energy. Thus the 
emition and absorption of phonon system by electric external field is essential for photo- 
induced effect. An external electric field leads to the catastrophe for the domino system. 
The paraelectric comes from quantum fluctuations between dominoes. This just the thermal 
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FIG. 2: Time dependence of electron order parameter. (3^1 = 0.01 and 0.1 aco phonon 

fluctuation induced quantum domino effect. The displacement of the dipole is suppressed 
because of the repulsion between the dipole moment. 

Next, we divide system A and B sub lattice. The Bloch-Wigner representation has to 
generate effective energy. First, the property of Liuville operator is given by. 



dt n 



{W,,g''[x,p]}, (9) 
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dx dp dp dx ' 
Wa = -Ito cos p + h'^. (11) 

The obtained effective energy is derived Moyal-Krammer formulation, as follows, 

^e//= (12) 
(T,a 

where a and A or B. Present representation of effective energy 

^e// = <ni> -2to cosp) <n^> . 

a 

< > -2to cosp) < > . (13) 




FIG. 3: Time dependence of spin order parameter. (3ei = O.Oland 0.1 no phonon 

where we neglect the fock term. Next we assume the system is Wigner crystal. The charge 
density is Akp and spin density is non magnetic. 

< K >= < >= (14) 

Thus, effective energy is given by 

1 - A2 

Eeff = U — 2tocosp (15) 

o 

Then we examine Landau expansion which obtained by free energy, 

^e// = ^e//|A=o + ^A + ^^A^ + ^BT J dxdp[f]nf+{l-f)]n{l-f)]. (16) 
The last part of the above equation is approximated as follows, 

J dpdx\hiY^ 
+ keT J dpdx[j + Y^]\A=o{^f\A=o^ 
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FIG. 4: Time dependence of spin order parameter. P^i = 0.01 and 0.1 aco phonon 



+ ksT J dpdx 
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In summary 2-nd order expansion is 

^2 = j(l - jp(eF)), 

where p(eF) is density of state of electron. The 3-rd order coefficient is 

_ kBTp{eF)U 

^."^ — : . 



The expression of 4-th order is given as 



kBTp{eF)U 



(17) 
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Thus, this system is asymmetric double well system. However, due to the low dimensional 
fluctuation there are no phase transition. 

In summary, we found a thermally excited paraelectric-to-ferroelectric crossover. We ex- 
amined time evolution of electron system with spin degree of freedom and phonon system. 
The quantum noise corresponds phonon excited by the electric external field. The imagin- 
ably random initial condition corresponds to emition and absorption of excited phonon by 
electric external field. The stronger randomness leads to metastable spin pierls state. The 
metastable state is unstable non-equilibrium state. This is ferroelectric state. We suggest 
such a non-equilibrium effect is necessary to photo induced effects. We also examine the 
Landau expansion of electron system at quarter-filled. 



* Electronic address: |it akurat@s6. dion.ne.jp 

^ K. Nasu: Relaxations of Excited States and Photo-Induced Phase Transitions: Proceedings of 
the 19th Taniguchi Symposium Kashikojima, Japan. Taniguchi International Symposium on the 
Theory of Condensed Matter 19. Springer Verlag. (1997). 

2 M. Takesada, T. Yagi, M. Itosh and S. Koshihara: J. Phys. Soc. Jpn. : 72 37 (2003). 

3 H. Seo and H. Fukuyama : J. Phys. Soc. Jpn. 66 1249 (1997). 

^ W. H. Press, S. A. Teukolsky, W. T. Vettering, B. P. Flannery: Numerical Recipes in C, Cam- 
bridge Univ. Pr. 2nd (2002). 



9 



